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ABSTRACT 

Many sulphur-bearing species have been detected in different astronomical 
environments and have allowed to derive important information about the chem¬ 
ical and physical composition of interstellar regions. In particular, these species 
have also been showed to trace and probe hot-core environment time evolution. 
Among the most prominent sulphur-bearing molecules, SO, sulphur monoxide 
radical, is one of the more ubiquitous and abundant, observed also in its isotopic 
substituted species such as 34 SO and S 18 0. Due to the importance of this simple 
diatomic system and to face the challenge of modern radioastronomical facilities, 
an extension to THz range of the rare isotopologues of sulphur monoxide has been 
performed. High-resolution rotational molecular spectroscopy has been employed 
to extend the available dataset of four isotopic species, SO, 34 SO, S 4 '0, and S 18 0 
up to the 1.5 THz region. The frequency coverage and the spectral resolution 
of our measurements allowed a better constraint of the molecular constants of 
the four species considered, focusing especially for the two oxygen substituted 
isotopologues. Our measurements were also employed in an isotopically invari¬ 
ant fit including all available pure rotational and ro-vibrational transitions for all 
SO isotopologues, thus enabling accurate predictions for rotational transitions at 
higher frequencies. Comparison with recent works performed on the same system 
are also provided, showing the quality of our experiment and the improvement 
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of the datasets for all the species here considered. Transition frequencies for this 
system can now be used with confidence by the astronomical community well 
into the THz spectral region. 

Subject headings: ISM: molecules — line: identification — molecular data - 
molecular processes — radio lines: ISM 


Introduction 


Molecules with sulphur account for about 10% of the species identified in the interstellar 
gas and circumstcllar envelopes, and among them the most prominent is sulphur monoxide 


(SO). Since its first detection by Gottlieb & Ball (1998) towards Orion A, many astronom¬ 


ical studies have been dedicated to this molecule, including those on the gas kinematics 


(Plambeck et al. 1982), molecular abundances (Blake et al. 1987), and spatial distribution 


(Sutton et al. 1995). Sulphur-bearing species are considered good probes of hot-core time 


evolution, since they are highly sensitive to physical and chemical variations of these regions 


(Viti et al. 2001 Hatched et al. 1998). Thus, they can be used as tools for investigating the 


chemistry and physical properties of complex star-forming regions (SFRs) located in dense 
molecular clouds. On the other hand, it is known that some molecules (SiO, H 2 CS, SO, S0 2 ) 


show increased abundances in regions affected by shocks (Bachillcr 1996) as a result of the 


action of outflows on the surrounding gas. The study of molecular lines from shocked areas 
provides valuable information about the chemical processes and the physical conditions of 
the shocked components. 


Recently Esplugues et al. (2013) dedicated an entire study to the sulphur oxide species 
in Orion KL, detecting 68 lines of SO, 34 SO, 33 SO and S 18 0 and providing upper limits 
for abundance of rarer isotopologues, such as S 17 0, 36 SO, and 34 S 18 0. Making use of the 
IRAM-30m antenna and analysing the SO non-LTE excitation they found that the most of 
the emission of SO arises from the High Velocity Plateau (column density N = (5 ± 1) x 
10 16 cm -2 ), and the Hot Core. These results confirm how SO (and S0 2 ) is a good tracer 
not only of shock-affected areas, but also of hot dense gas. The analysis of the abundance 
ratio of the SO isotopologues compared to that of S0 2 for the different component of the 
observed region, showed that the latter is five time more abundant than SO in the High 
Velocity Plateau and the Hot Core. This work clearly elucidated the importance of isotopic 
analysis in the interstellar medium to derive information both on the abundance of optically 
thick parent species, and on the evolution of the isotopic ratio in different astronomical 
environments. 
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More recently, a new extensive laboratory spectroscopic characterisation of sulphur 
monoxide isotopologues has been published by Martin-Drumel et al. (2015, hereafter MD). 
In this work the authors provided new high-frequency measurements for the main isotopic 
species and its sulphur isotopic substituted ( 34 SO and 33 SO) using continuous-wave terahertz 
photomixing based on a frequency comb, with a frequency accuracy ranging from 10 to 800 
kHz. Furthermore, the authors combined their measurements with all available pure rota¬ 
tional and ro-vibrational transitions of different SO isotopologues in an isotopically-invariant 
fit. 


In the work here presented, we report the extension to the THz-range of the available 
dataset for sulphur monoxide and the minor isotopologues 34 SO (4.2%), S 18 0 (0.2%) and 
S 17 0 (0.04%) through laboratory absorption spectroscopy, in order to improve the available 
spectroscopic parameters and determine centrifugal distortion terms. The determination of 
the latter has also supported by high-level quantum-chemical calculations. Finally, a detailed 
comparison with previous available spectroscopic data is supplied in the discussion section, 
together with the results from an isotopic invariant fit involving all available rotational and 
ro-vibrational frequencies for all isotopologues of SO. 


2. Experiment 


All the measurements have been carried out with a frequency-modulated, computer- 
controlled spectrometer (65 GHz - 1.6 THz) (Cazzoli & Puzzarini 2006 Puzzarini et al. 
2012 Cazzoli & Puzzarini 2013), with the 200 GHz - 1.5 THz range actually considered. 


Frequency multipliers driven by Gunn diode oscillators, phase-locked to a rubidium frequency 
standard, have been employed as the millimeter and sub-millimeter-wave sources. The fre¬ 
quency modulation was obtained by sine-wave modulating the 72 MHz local oscillator of 
the synchronisation loop at 16.66 kHz, with the modulation depth varied from 400 kHz 
to 2.4 MHz according to the transition frequency under consideration. Liquid-He-cooled 
InSb and Schottky diode detectors were used and their output processed by means of a 
lock-in amplifier tuned to twice the modulation frequency (i.e., second harmonic detection 
was performed). Previous unrelated experiments performed in the cell with sulphur-bearing 
chemicals allowed us to achieve very good signals of sulphur monoxide only with the flow of 
~ 30 mTorr of O 2 and a DC discharge of 12 mA, at room temperature. The measurements 
of the 34 S isotopic species have been carried out in natural abundance (4.22%), while an iso- 
topically enriched oxygen sample, a mixture of 10% 17 02 and 90% 1S 02, has been employed 
for the experiment of the two oxygen minor isotopologues. A few-Gauss magnetic held, par- 
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allel to the radiation source, has been used to disentangle the free-radical signal from those 
of closed-shell species whose transitions contaminated the spectral window. Subtraction of 
the on-field spectrum to the off-one allowed a perfect zeroing of the effect of the undesired 
non-magnetic species (see Fig. [Tj). 


3. Computations 


As mentioned in the Introduction, to guide the analysis and the experimental deter¬ 
mination of the sextic centrifugal distortion constant, quantum-chemical calculations at the 
coupled-cluster (CC) singles and doubles (CCSD) approach augmented by a perturbative 


treatment of triple excitations (CCSD(T); Raghavachari et ah 1989) were performed with 
the cc-pCVQZ basis set (Dunning |1989 Woon & Dunning 1995 |Peterson et al. 2002). All 
electrons were included in the correlation treatment. 


For each isotopic species considered (SO, 34 SO, S 17 0, S 18 0), the harmonic force held was 


obtained using analytic second derivatives (Gauss & Stanton 1997), whereas the cubic force 


held was determined in a normal-coordinate representation via numerical differentiation of 


the harmonic force constants (Schneider & Thiel 1989 Stanton et ah 1998). The quartic 


and sextic centrifugal-distortion constants were then determined by means of vibrational 


perturbation theory (Mills 1972 Puzzarini et ah 2012). 


All computations were performed with the CFOUR program package (2012). 


4. Analysis 

The total angular momentum J for a diatomic molecule in the 3 E state is given by the 
coupling between the molecular angular momentum N and the electronic spin momentum 

S. The effective Hamiltonian operator may be expressed as the sum of three terms: 

H = H rot + H S -n + H S -s 

with 

H rot = N 2 (B- DN 2 + FN 4 ) 

being the rotational Hamiltonian. For each N value (being N the associated quantum num¬ 
ber to the angular momentum operator), there are three different energy levels corresponding 
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to the three possible space orientation of the S vector (5=0,±1). 

Hs~n = N ■ S(7b + 7dN 2 ) 


H S -s = 1(\b + \dN 2 )(3S 2 z -S 2 ) 

are Hamiltonians describing the electron spin-rotation and spin-spin interactions, respec¬ 
tively. Bj'Jbj and A b are the rotational, the electron spin-rotation and the spin-spin interac¬ 
tion constants, and D,H, 7 ^, and A d represent their centrifugal distortion correction. The 
derived parameters of the Hamiltonian for each isotopic species are summarised in Table [l](2j 

The effective Hamiltonian of the 1( 0 species also includes the nuclear hyperfine inter¬ 
actions with the overall rotation, due to the 17 0 nuclear spin (/ = 5/2): 

1 (‘U 2 —I 2 1 

H hfs = b F I ■ S + c(I z S z - -I ■ S) + eQq ^ * _ ^ + CVI • N 

where bp and c are the isotropic (Fermi contact interaction) and anisotropic parts of the 
electron spin-nuclear spin coupling, eQq is the nuclear quadrupole coupling constant, and 
Ci is the nuclear spin-rotation interaction constant. 


Sulphur monoxide has been the subject of several laboratory studies. The Cologne 
Database for Molecular Spectroscopy (Muller et al. 2005) has been used as a review of the 
previous works, and to obtain a starting dataset to flag transitions to remeasure with higher 
accuracy or to derive higher frequency predictions. The first spectroscopic study for the 
fundamental isotopic species date back to the late 70s (Clark & De Lucia 1976), followed 
in the years by many other works, including some THz measurements (Cazzoli et al. 1994). 
Also for the 34 SO previous laboratory data cover quite a large range of the mm/sub-mm 
wave spectrum, with some measurements taken at 1.05 THz by (Klaus et al.||1996), allowing 
a total of 42 rotational transition of the X 3 T,~ electronic ground state catalogued. Somewhat 
more critical have been in the past years the measurements for the minor oxygen isotopic 
substituted species, such as S 18 0 and S 17 0. Both species have been the subject of the same 
work of the 34 S species, but while for the former the measurements have been carried out 
up to the THz region (1.03 THz), the rarer 1( 0 species has been studied up to 620 GHz. In 
all the cases mentioned above, the accuracy of the measurements spans from 50 kHz, in the 
lower frequency range, up to to 100 kHz in the THz region, values that our measurements 
can easily improve. 


As briefly mentioned in the introduction, during the conclusion of the present study, 


a new spectroscopic analysis on the SO isotopologue system has been published (Martin- 


Drurncl et al. 2015). The focus of this work has been the extension of the absorption 
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spectrum of SO, 34 SO, and 33 SO up to 2.5 THz, allowing to perform a global isotopically- 
invariant fit.In the present work, for the main isotopic species, two rotational transitions in 
the 340 GHz region have been remeasured and 10 new transitions have been detected in the 
1.07-1.41 THz frequency range. The regions 200-500 GHz and 1.10-1.40 THz have been 
studied for detecting 23 and 24 rotational transitions of the 34 SO and S 18 0 isotopic species, 
respectively. Finally, the rarest isotopologue S 4 '0 has been observed through 15 rotational 
transitions in the 1.10-1.50 THz region. The accuracy of our measurements is 20-30 kHz 
in the millimeter-/submillimeter-wave range and 50 kHz in the THz region. The measured 
frequencies will be made available for the spectroscopic community through the CDMS on 
line catalogue. 


5. Discussion 

The results of these new measurements are summarised in Tables [l] and [2j for SO and 
34 SO (rms = 0.62 and 0.86), and S 17 0 and S 18 0 (rms = 0.74 and 0.54), respectively, and 
compared with the data obtained by MD. From this comparison we observe that our labo¬ 
ratory data for SO and 34 SO are in a very good agreement with the corresponding constants 
derived from an effective-fit for the same species (see Table 2 of MD). The higher frequency 
measurements performed in the previous work for the main species allowed the authors to 
constrain the second-order centrifugal distortion correction of the spin-spin interaction con¬ 
stant, A 2 , undetermined in our analysis. All other terms show an agreement of better than 
1 %, except for the sextic centrifugal distortion parameter, H: in this case our parame¬ 
ters are larger than those derived previously by 11 and 19%, respectively for SO and 34 SO. 
To inspect this large deviation, we make use of our computed values and of an empirical, 
well-tested scaling procedure: 


TTISO _ TT 

n seal ~ n ‘ 


calc 


X 


( Tjvnain / Tjmain\ 
\-H-exp / n calc ) 


where the superscripts iso and main refer to a specific isotopic species and the main iso¬ 
topologue, respectively; seal, exp, and calc denote the scaled, experimental, and quantum- 
chemically calculated values for H , respectively. Using the value of -6.974(22) x 1C) -9 MHz for 
SO, we obtain a scaled value of -6.56xlCD 9 MHz for 34 SO, which is roughly in between our 
and MD values and well within the confidence range of our datum. This outcome suggests 
that the fit of MD might take some benefit from the inclusion of the present measurements. 
In Table [2] the results for the two oxygen-substituted species are collected. Here, the liter¬ 
ature values reported for the comparison are those predicted from the isotopically-invariant 
fit available in MD. As for the previous table, the agreement for all the parameters is bet- 
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ter than few percent. The only relevant deviation is noted for the H constant of the S 17 0 
species: -6.31(28) xlCD 9 MHz vs -5.74xl0“ 9 MHz. Once again, we resort to the scaled com¬ 
puted value (-6.18xl0 -9 MHz) to a deeper understanding of such a discrepancy. The latter 
value well agrees, within the given uncertainty, with our datum, while it suggests that the 
H constant from the isotopically-invariant fit is underestimated. Therefore, this fit would 
largely benefit by the inclusion of our recorded frequencies. Since the Dunham expansion 
model MD used to determine the predicted parameter does not include hyperfine interac¬ 
tions, in Table [2] the corresponding parameters are compared with the results from Klaus 


et ah (1996). A very good agreement is observed. The only comment concerns the 11 0 


spin-rotation constant, which is badly determined by the fit: we decided to keep it fixed at 
the quantum-chemical value, computed at a level of theory known to provide quantitative 
agreement with experiment (Puzzarini et al. 2010b). For S 18 0, we notice that the scaled 
value of H is -5.53 xl0~ 9 MHz and both our and MD values seem to be underestimated. On 
the other hand, our datum is affected by a large uncertainty and the scaled value is within 
the upper limit. 


It is also worth to compare now, for the oxygen isotopologues of sulphur monoxide, all 
the three available datasets, specifically ours, the MD one, and the catalogue on the CDMS 
website ( www.cdms.de ). In Figure [2] this comparison is made for the predicted transitions 
in the frequency range 60-2000 GHz; in particular, it is illustrated the difference of the 
strongest set of transitions ( AN,AJ = +1, with J = N, N ± 1), the more relevant for 
the astrochemistry community, among our dataset and the two previous ones. The need of 
new measurements is particularly evident when comparing our new predictions with those 
available on the CDMS catalogue, which, in the high-frequency range (> 1 THz) differ by 
more than 1 MHz from ours, for the 17 0 isotopic substituted species. Less significant but 
still sizeable is the comparison with the Dunham isotopically-invariant results. First, it is 
important pointing out the overall quality of these predictions, which make use of a very large 
dataset of measurements of different isotopic species and several vibrational excited states. 
Second, it must be noted that also in this case, although less prominent than in the previous 
one, the deviation at high-frequency might be worth to consider, and very likely related to 
the discrepancy observed for the sextic centrifugal distortion term. The two datasets show 
a considerable better agreement for the other oxygen species, while a larger difference can 
be noted comparing our predictions with the older data. 


To conduct a more exhaustive comparison with the previous data and to inspect the 
impact of our new measurements, a multi-isotopologue Dunham analysis of the full dataset 










has been performed, using the known reduced mass dependences given by: 


Yi, m = Ui 


for the rotational and vibrational parameters; 


X, m = U, 


with X = A, bp, c, eQq-, and 



-\-m 



X\ m = Ui 


X 

l.m 



with A" = 7,(7/. The A\ m and Aare unitless coefficients which account for the effect of 
Born-Oppenheimer breakdown. For a clearer relation of the Dunham’s coefficients with the 
molecular energy levels and the rotational-vibrational parameters see Townes & Schawlow 


(1975) equations 1-35 and 1-36. The SPFIT/SPCAT software (Pickett 1991) has been 


used for this analysis and the corresponding parameters are shown in Table [3j The SPFIT 
parameters have been converted to the isotopically invariant terms U^ m through the m 
and X i m relations aforementioned, for a better comparison with the previous published 
data; more specifically the parameters derived from the present work are set side by side 
with the ro-vibrational isotopically constants of MD, along with the relative reduced rms 
errors. Some high-order parameters derived in MD have been left out from our analysis 
since their inclusion in the global analysis did not contribute to the quality of the final 
fit. The comparison shows a general very good agreement for all the parameters analysed. 
Nevertheless, there are some discrepancies that required a deeper look into. The vibrational 
corrections to the rotational constants, [Ayi, tAu, and U$a, show a difference from a factor few 
to more than 10. Surprisingly, both the former two parameters are in very good accordance 


(in la deviation) with the isotopically invariant analysis performed by Klaus et ah (1996) 


to which MD compared their work (see Table 4 of MD). Unfortunately this comparison can 
not be addressed for the third parameter U- h \, which have been determined for the first 
time by MD. A similar situation can be seen for the higher order vibrational corrections to 
the A constant, namely the U£ Q , U§ 0 and U$ 0 parameters. The former constant is the one 
showing the largest deviation, including a change in sign, from -0.1057 to 0.6547 MHz-amu 2 . 
Again, as for the previous cases, a more accurate analysis reveals a good accordance (in la 


deviation) with the value derived by Klaus et al. (1996). Similarly, the other two parameters, 


Kg 0 and f/g 0 , have been determined for the first time by MD and can not be compared to 
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a third analysis: here the difference is less critical among the two studies but there is still a 
factor ~4 and ~ 10, respectively. 


Crucial parameters of an isotopically invariant analysis are the Born-Oppenheimer 
Breakdown constants. In Table [4] the BOB constants derived from our work are compared 
with those from MD. The parameters for the correction to the rotational constants, Aq X 
and Aq X , show a very good agreement with the latter slightly (~ x2) improved with our 
global analysis. Similar good accordance is found for the sulphur correction to the vibra¬ 
tional frequency, Af 0 , while a large discrepancy, including a change in sign, is shown for 
the oxygen correction, Af 0 . A future study including more ro-vibrational and vibrational 
transitions of oxygen substituted species (currently there is only the first excited vibrational 
band of S 18 0 in the analysis) might help to elucidate this discrepancy and further constrain 
this important parameter. To compare the BOB coefficients of the fundamental A constants, 
Aq’q and Aq’q, the equation 



K 

+ Xfa + X 0 B 0 ) 



has been used. Here A is the atom considered, M^ is the atomic mass of the main isotopic 
species of the A atom, m e is the electron mass, U is the Dunham constant to which the 
BOB coefficient is related, and the X’s are the parameter provided from MD analysis. The 
oxygen BOB coefficient is in very good agreement with that derived in the previous work, and 


slightly improved (~ x2); the sulphur one, Aq’q , is dissimilar and, once again, opposite i 
sign. As seen previously for the discussion about the I/’s, a comparison of our analysis with 
that of Klaus et al. (1996) reveals a better accordance for this parameter, with respect to its 
sign (now both positive) and its absolute values, although different to few sigma deviation 
(0.886 and 0.221, our and Klaus et al. 1996, respectively). 


m 


The present laboratory work extends the spectroscopic study of a very simple and im¬ 
portant diatomic radical well into the THz region, including some minor isotopic species. 
For such an important molecular system, widely distributed and very abundant in many 
different interstellar regions, the accuracy of the dataset of less abundant isotopologues is 
also crucial. This is particularly true in an era in which high-sensitivity, high-spectral and 
spatial resolution radioastronomical facilities are becoming available to the scientific com¬ 
munity. The ALMA interferometer has become the reference and most relevant observatory 
for the mm- and sub-nnn-wave region; its unprecedented properties will allow to observe SO 
and its isotopologues in regions previously undetected. Above the THz threshold the instru¬ 
ments of the airborne observatory SOFIA will be also able, in the near future, to possibly 
extend the identification of this molecular system in the interstellar environments at higher 
frequency. 
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Fig. 1.— The N,J = 9,10 8,9 rotational transition of 34 SO. The black-dashed shows the 

first acquired spectrum, in which a contaminant species transition is on the left edge of the 
34 SO rotational transition. The effect of applying a magnetic field is to cancel the radical 
transition, isolating the non-magnetic contaminant transition (in red), whose contribution 
can then be easily subtracted to derive a precise information on the radical of interest (in 
blue). Owing to the modulation and detection scheme employed, the instrumental lineshape 
is approximately the second derivative of a Lorentzian. No averaging time was required for 
obtaining a good S/N. 
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Fig. 2.— Comparison of the predictions (< 2 THz) of the AN, A J = +1 (with J — N, Ail) 
transitions, for the two oxygen isotopic substituted species. For each N the three J levels 
are shown. “New” predictions are based on our data set. 




















Table 1. 


Spectroscopic constants' 1 of SO and 34 SO (in MHz) 



SO 

34 SO 

This work 

MD 6 

This work 

MD b 

B 

21523.55568(26) 

21523.555878(79) 

21102.73094(43) 

21102.73228(31) 

10 3 D 

33.91431(60) 

33.915261(85) 

32.59862(90) 

32.60073(55) 

10 9 H 

-7.73(37) 

-6.974(22) 

-7.22(49) 

-6.06(31) 

As 

158254.3860(80) 

158254.392(09) 

158249.715(12) 

158249.807(24) 

Ad 

0.30662(12) 

0.306259(72) 

0.30030(38) 

0.30063(08) 

10 6 Ad 

- 

0.478(47) 

- 

- 

7 B 

-168.3045(22) 

-168.3043(15) 

-164.9863(33) 

-164.9948(36) 

10 3 7d 

-0.5239(46) 

-0.52545(97) 

-0.5181(61) 

-0.5116(27) 

a 2 

0.62 

0.88 

0.86 

0.85 


Uncertainties (lcr) are in units of the last significant digit. 


b Martin-Drumel et al. (2015). 
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Table 2. Spectroscopic constants' 1 of sulfur monoxide oxygen minor isotopologues (in 

MHz) 



O 

00 

GO 

s 17 

O 

This work 

MD b 

This work 

MD b 

B 

19929.27860(42) 

19929.27880 

20677.80826(76) 

20677.80911 

10 3 D 

29.06925(96) 

29.0683 

31.29739(82) 

31.2977 

10 9 H 

-5.10(60) 

-5.14 

-6.31(28) 

-5.74 

As 

158236.114(40) 

158236.063 

158244.81(26) 

158244.75 


0.283118(39) 

0.283426 

0.29555(26) 

0.29413 

7 B 

-155.8017(61) 

-155.7908 

-161.660(24) 

-161.666 

10 3 7 d 

-0.4366(78) 

-0.4510 

-0.4889(85) 

-0.4855 

b F 

- 

- 

-39.784(17) 

-39.783(ll) c 

c 

- 

- 

95.359(30) 

95.340(48) c 

eQq 

- 

- 

-3.737(78) 

-3.614(60) c 

10 3 cv 

- 

- 

-7.57 d 

-4.4(15) c 

a 2 

0.74 

_e 

0.54 

- 


Uncertainties (lcr) are in units of the last significant digit. 


b Martin-Drumel et al. (2015). 
c Klaus et al. (1996). 

d Kept fixed at the computed value. See text. 


e 


Spectroscopic parameters are computed values derived from the isotopically invariant 
fit. 










Table 3. Isotopically invariant parameters 0 for the SO radical 


Parameter 

This work 

MD 

Units 

t/i,o c 

3757.50375(287) 

3757.49542(195) 

cm -1 -amu 1 / 2 

U v 0 C 

-68.32282(682) 

-68.32167(622) 

cm 1 -amii 

U 3 , 0 C 

0.47701(898) 

0.47556(818) 

enr'-amu 3 / 2 

C/ 4 ,o c 

-0.04119(509) 

-0.04041(463) 

cm -1 -amu 2 

C/ 5 ,o c 

0.00519(104) 

0.005031(934) 

cm -1 -amu 5 / 2 

U 0 ,i 

230416.7278(172) 

230387.59070(231) 

MHz-amu 

c/ M 

-6001.5217(148) 

-6001.6578(164) 

MHz-amu 3 / 2 

U 2 ,i 

25.4579(211) 

25.7585(267) 

MHz-amu 2 

C/ 3,1 

-0.5585(120) 

-0.8254(175) 

MHz-amu 5 / 2 

C/ 4 , 1 

-0.12215(270) 

-0.01197(508) 

MHz-amu 3 

c/5,1 

-0.001253(218) 

-0.022075(633) 

MHz-amu 7 / 2 

c/0,2 

-3.8542930(244) 

-3.8543097(229) 

MHz-amu 2 

c/1,2 

-0.004668(126) 

-0.004547(114) 

MHz-amu 5 / 2 

c/2,2 

-0.0010726(628) 

-0.0011636(577) 

MHz-amu 3 

c/0,3 

-7.8512(835)E-06 

-7.8489(791)E-06 

MHz-amu 3 

c/0,4 

-1.023(156)E-09 

-1.040(146)E-09 

MHz-amu 4 

c/0,0 

-1787.3375(287) 

-1787.4277(276) 

MHz-amu 

c/7,0 

-45.818(188) 

-45.508(176) 

MHz-amu 3 / 2 

Ulo 

2.0778(143) 

1.871(133) 

MHz-amu 2 

TTl 

^3,0 

0.0721(227) 

0.0986(216) 

MHz-amu 5 / 2 

C/ 0 7 ,1 

-0.060730(198) 

-0.060460(184) 

MHz-amu 2 

c/7,1 

0.00454(117) 

0.00439(108) 

MHz-amu 5 / 2 

77 A 

( - 7 o,o 

157785.078(303) 

157795.4203(201) 

MHz 

U x 

ty i,o 

2978.880(153) 

2979.597(136) 

MHz-amu 1 / 2 

U x 

u 2fl 

114.524(188) 

112.732(159) 

MHz-amu 

TJ X 

u 3,0 

11.905(107) 

13.6225(953) 

MHz-amu 3 / 2 

IJ X 

^4,0 

0.6547(290) 

-0.1057(327) 

MHz-amu 2 

7V A 
^5,0 

0.05997(282) 

0.22948(606) 

MHz-amu 5 / 2 

TJ X 

^ 6,0 

-0.001964(177) 

-0.018122(580) 

MHz-amu 3 

U X 

u 0 ,l 

3.24178(138) 

3.24222(125) 

MHz-amu 

TT X 

u l,\ 

0.14168(874) 

0.14226(791) 

MHz-amu 3 / 2 
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Table 3—Continued 


Parameter 

This work 

MD 

Units 

U x 

u 2,l 

0.03660(511) 

0.03545(466) 

MHz-amu 2 

U x 

U 0,2 

5.626(387)E-05 

5.374(357)E-05 

MHz-amu 2 

IJ bf 

u 0,0 

-39.7904(214) 

-39.7830(110) 

MHz 

TJ C 

^0,0 

95.3590(315) 

95.3400(480) 

MHz 

1.5 X U e Qq 

-3.5687(859) 

-3.6140(600) 

MHz 

o 

X 

■ft 

1 

o 

-0.00757 

-0.00440(150) 

MHz 

a 2 

0.7491 

0.8997 



a Uncertainties (la) are in units of the last significant digit. 


Martin-Drumel et al. (2015). 


c The units for these constants have been corrected: in MD these parameters were 
missing the amu dependencies. 


Table 4. Born-Oppenhcimer Breakdown coefficients' 1 derived from the fit 


Parameter 

This work 

MD b 

Klaus et al. c 

A s 

*-*0,1 

-1.9737(16) 

-1.9725(14) 

-1.9772(58) 

A° 

^0,1 

-2.7015(20) 

-2.7175(57) 

-2.7247(34) 

A 5 

^1,0 

0.409(35) 

0.395(32) 


A° 

^1,0 

-0.2664(89) 

0.1338(40) 


A A.5 

^0,0 

0.866(55) 

-0.110(85) 

0.211(89) 

A a ’° 

^0,0 

1.494(46) 

1.427(87) 

1.330(48) 


“ Uncertainties (la) are in units of the last significant digit. 
b Martin-Drumel et al. (2015). 


Klaus et al. (1996). 
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